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ABSTRACT 


1. Introduction 


In the last few years, theory and application of PIDEs play an important role in the various fields of science and 
technology. It is very important to know various methods to solve such partial integro differential equations. It is 
found that, the double Laplace transform [1] and Double Elzaki Transform [2] used to solve PIDEs. Using DLT and 
DET, linear PIDEs convert into an algebraic equation then after factorizing it and applying their inverses an exact 
solution of the problems was achieved. Recently, Modified Differential transform method [11] and two 


dimensional DTM [10] have used to solve linear PIDEs. 


Laplace Differential Transform Method is the combination of the Laplace Transform method and one-dimensional 
Differential Transform Method. Earlier, the LDTM has been used to explore the exact solutions of Linear and 
Nonlinear Klein Gordon equations [5]. Similarly, the same method has used for nonlinear PDEs with boundary 
conditions [6] and Wave equations [7]. According to the latest research, LDTM has used to solve nonlinear non 


homogeneous Partial Differential Equations [8]. 


In this article, first convert linear Partial Integro Differential equation with convolution kernel into Ordinary 
differential equation, then solved this ODE by Differential transform method and then an exact or series solution is 


achieved by applying inverse Laplace transform. 

2. Preliminaries 

2.1. Partial Integro Differential Equation 

The general form of PIDE with convolution kernel is, 

(dtu ~ atu . : atu 

2 aj (= + Za bj (Fr) + cu(x,t) + 2 a, | K,(t — y) (= dy + f(x,t) =0 (1) 
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Where, a; ,b;,c and d; are constants or the functions of x alone and f (x,t) , K;(t — y) are known functions. 


2.2. Laplace Transform 


The Laplace transform of a function f(t) is define by, 


f(s) = Lf} = e“tp(tdt ; t20 
0 


Where, s is real or complex. 
The Laplace transform of first and second order partial derivatives are as follows: 


Ou 


i) Lf} = Uu,(%,S) 


. (I : 
ii) Lis = Uyy (X, 5) 


Ox? 
Ou _ 
iit) L ea = su(x,s) — u(x,0) 
: d*u a 
iv) L {Fat = s*iu(x,s) — su(x,0) — u;(x, 0) 


Using Laplace transform convolution theorem becomes, 
7 t 
Lg(s).f(S)} =fO «9 = fmat —y)dy 
0 


2.3. Double Laplace Transform (DLT) [1] 


The Double Laplace transform of a function f(x, t) in the positive quadrant of the xt-plane is defined by, 
foe eL,14Afa) = | ev | eS*f(x,t)dxdt ; x,t>0 
0 0 


Where, p,s are real or complex. 


The double Laplace transform of first and second order partial derivatives are as follows: 
of = - 
i) Lele |} = pfs) - Fs) 


(a) _ S 
ii) Lyly ny = sf(p,s) — F(p,0) 


i a OD ioe paick 
iit) Leb {rat =P F(@,s) — PFO,s) - £5) 
, 8 ee eae 
iv) LyLy {Fart = s*F(p,s) — sf, 0) — fe, 0) 
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2.4. Double Elzaki Transform (DET) [2] 


The Double Elzaki transform of a function f(x, t) in the positive quadrant of the xt-plane is defined by, 


T(p,s) = E,E,{f (x, t)} = ps | [ve “G a ) F(x, t)dxdt ; x,t>0 


Where, p,S are real or complex. 


The Double Elzaki transform of first and second order partial derivatives are as follows: 


0 1 
i) ExEe {5} ==10,5) - pT(0,5) 


0 1 
ii) EF, ea = si @, s) — sT(p,0) 


i 1 
iii) E,E;, [os | = pe s) —T(0,s) — pT,(0,s) 


a 1 
i) Exbi {Soa} = S3T.8) 7.0) ~ 5M.) 


2.5. Differential Transform Method (DTM) [3,4] 


Let the function u(x) is arbitrary analytic and continuously differentiable function of x defined on D = [x,0] ¢ 


R and x9 € D then, 


1 ( a* 
‘ {x=Xo} 
. u(x) = > U(k)\(x—x)* ; k 20 


Note: If x9 =0 then ”~ u(x) = > U(k)x* ; k>0 
k=0 


We give the proofs of properties of one dimensional DTM in the following theorem: 
Theorem 1: If U(k) is differential transform of u(x) then 


k=m 


a) If u(x) = x™ then U(k) = 6(k —m) =| mau 


k 
a 
b) If ux) = e™ then U(k) = a 


a® kn 
c) If u(x) = sinax then U(k) = —sin(> 3 ) 
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d) If u(x) = cosax then U(k) = 70 £5 (=) 


e) If u(x) = ae then U(k) = (kK +1)(K+2)...(k +m)Y(k +m) 


k 
f)If u(x) = f()g() then Ud) = > F(m)G(k —m) 
m=0 


Proof: a) Let, u(x) = x™ 


ak* 
x=0 
= {m(m - 1)...(k —-m+1)x™*},_9 
If k=m -U(k)=1 
lf kem UG) =0 
1 ; 
& UG=6h a= f 
b) The proof is obvious. 


c) Let, u(x) = sinax 


ak 
“ k!U(kK) = ta sinax| 
x=0 


=fatsin(ax+ 2) -ue) = Sin“) 
Similarly we can prove (d). 
e) Let, u(x) = 2 
= Uda = [re Say 
=(k+m)!Y(k+m) » U(k) = (kK +1).(k + 2)..(k + m)Y(k +m) 
f) Let, u(x) = f)g@) 
ah 
“ KLU(k) = [5 renga} 
-| : aes stio= y F(m)G(k —m) 
a on oar os a 


We summarize the properties of one dimensional DTM in the following table: 
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Table 1. Properties of one dimensional DTM 


Original Function u(x) Transformed function U(k) 
a hei) 
: 0; k#0 
1; k=m 
m ) 
" ti ; k#m 
eX ae 
k! 
sinax i 
wz) 
cosax at ud 
mes) 
0 
a (k + 1Y(k +1) 
Ox 
amy ig ee ete Nee (ic om) (ic 
yee +m) 
k 
f(x) 9g) >. F(m)G(k — m) 
m=0 


3. Applications to Linear PIDE 
3.1. Double Laplace Transform Method 


Applying double Laplace transform to (1), 


m i-1 ai n i-1 ak 
pa a; |p'u(p,s) — Nas Lie (FH0 9) + 3 b; |s'u(p,s) — > ee (Fane 0) 
i= j= 


i=0 k=0 
Tr i-1 ai 
+eHp,s) +) di k(s)| > pit Ly (Fu 9] +f@,s)=0 
i=0 j=0 


This is an algebraic equation in u(p,s). After factorizing it and applying inverse DLT, an exact solution of (1) 


will be achieved. 
3.2. Double Elzaki Transform Method 


Applying double Elzaki transform to (1), 
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i-1 
a ok 
>, ai| 57,5) - Dor, Zr.) +a 218) 9's¥i8,(F7000)) 


i=0 k=0 


i-1 


Tr * 
= set d/ 
+ceT(p,s) + >» d, K(s) >, p2-i+i F, (Fre. 9) irc 
i=0 


J=0 
This is an algebraic equation in T(p,s). After factorizing it and applying inverse DET, an exact solution of (1) 
will be achieved. 
3.3. Laplace Differential Transform Method 
Applying Laplace transform w. r.t t to (1), 


>. a ee a+ oe “} retin) +) avafaen -$5] +-176501=0 


i=0 


SS atu SS atu . aiy 
* Qtegat* 2, Paget elas) + K(s) 5 
l= i= 


i=0 


}+ Fas) =0 (2) 


Equation (2) is an ODE. This ODE will be solved by one dimensional differential transform method and then by 


inverse LT, an exact solution of linear PIDE will be achieved. 
4. Results and Discussion 


Example 1: Consider PIDE 
t 
Ure = Uy + 2{ (t —a)u(x,a)da — 2e* (3) 
0 


With initial conditions u(x,0) = e*, uz(x,0) = 0 and boundary condition u(0,t) = cost. 
Double Laplace Transform [1]: Applying DLT, we get 


( a Vu a #8 : 2 Ss 
ae iS) a s(p-1) p-1 
s 


eC EEN CEED 


(4) 


Applying inverse DLT, 
u(x,t) = e*cost 
Double Elzaki Transform [2]: Applying DET, we get 


1 il ps? 2p’s p 
- 2s*- 5) Ss) =—— - ——____ - — 
(5+ o 7 (p, 5) s?+1 (1-p)(s?+1) 1-p 
22 
p’s 
7s) =>—+___—_ 
PS) = Goes) 
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Applying inverse DET, 


- u(x,t) = e*cost 
Laplace-Differential Transform method: 


Applying LT with respect to t on (3) and using u(x,0) = e*, u,(x,0) = 0 


OS) 4 (5 — st) ues) = e*(=-s) (6) 


Now applying DTM 


De +1,s) =~ {(5 =?) a +(- ) 7 
i ~ (k+1) s? ee 7) 
Now applying the Laplace transform to boundary condition 
U(0,s) = 2 U(0,s) = 2 
he) aed eS) aed 
Put k = 0,1, 2,3,.... in equation (7), 
fio. Fiojs— Ifk=1; U@ ==/ -_) 
— ae oe SPS DNs? £1 
fk=2; 0G,s)=-(= es) TE92 — d 
I = 2s U(3,s)=5(sz5) I = 3% 04,8) =F (az5) and so on... 
ae : ee 3 
Now, U(x, -) 0 k, = Hl _ 
ow (x, Ss) 2, (k, s)(x) aa et a hat 


Now applying Inverse LT, 
-. u(x,t) = cost e~ (Also shown in graphical presentation: see Fig. 1) 


Example 2: Consider PIDE 


t 
Up — Uy tut ) e!-Du(x, a)da = (x? + 1)e'-2 (8) 
0 


With initial conditions u(x,0) = x?, u,(x,0) = 1 and boundary conditions u(0,t) = t, u,(0,t) = 0 


Double Laplace Transform [1]: Applying DLT, we get 


1 2p 2. 2 P 
[sop tt GUO) = Ey pet st 


2 1 
LUGS) = 535 + (9) 


s ps? 
Applying inverse DLT, 
* Ut) = x* +e 


Double Elzaki Transform [2]: Applying DET, we get 
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1 1 s es oe 
~-S+1+——)r ,S)= (2p7 +1 — 2p’s* + 2sp* — s3 
(;-sati+7=,)P@9)= @p (=) p?s? + 2sp*—s 
aT (ps) = 2p" s* +p7s* (10) 


Applying inverse DET, 
» U(%,t) =x? +t 
Laplace-Differential Transform method: 


Applying LT with respect to t on (8) and using u(x, 0) = x?, u,(x,0) =1 


d?U(x,s) s? 2 (x7 +1) 
. U ==—=——__— =" 11 
ax? (5) us) s s-1 oo 
Applying DTM, 
U(k +2,s) = * Slates) ace 2)(—) 5) +280) 12 
Or (et + I I\e=1) c—/ s=4 5 -) 
Now applying the LT to boundary conditions we get, 
1 dU(0,s) 
UO,s)=>3 and ———=0 
Ss dx 
; = 1 = 
Applying DTM, “« U(0,s) = 52 and U(1,s)=0 
Put k = 0,1,2,3,... in equation (12), 
= 1 as 
If k=0 . U(2,s)=— ; If k=1 - U(3,s) =0 
fk=2 #0(4,5)=0 : Ifk=3 « U(5,s)=0  andsoon... 
i ta 
Now, U(x,s) = > U(k,s)(x)* = [+ —x? 
a s*  s 
-. Applying inverse LT 


“ u(x,t) =x?+t (Also shown in graphical presentation: see Fig. 2) 


Example 3: Consider PIDE 


t 
Ut + Utet — sinh(t — @) uyy,(x, ada = 0 (13) 
0 


With u(x,0) = 0, u;z(x%,0) =x, uge(x,0) = 0 and u(0,t) = 0, u,(0,t) = sint, u,,(0,t) = 0. 
Double Laplace Transform [1]: Applying DLT, we get 


3 
,  P s p 

_ U se ees 

(s ne s2 — ;) (p,s) p2 (s*-—1)(s74+1) 
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 U(p,s) = : 14 
UO) =p (14) 
Applying inverse DLT, 


. u(x,t) = xsint 


Double Elzaki Transform [2]: Applying DET, we get 


C 1 | )r — nes? 
s s* 1=p? eee (1 — p?)(s? +1) 
oj” 15 
TS) =3o5 (15) 
Applying inverse DET, 
. u(x,t) = xsint 
Laplace-Differential Transform method: 
Applying LT with respect to t on (7) and using u(x,0) = 0, u,(x,0) = x, uUpge(x, 0) = 0 
d3 
Es 7a UG s) —(s?-1)(s +. s3)U(x,s) + sx(s? —1) =0 (16) 
Now applying DTM 
= 1 = 
U(k + 3,s) = ——_.—__———_{(s° — s Uk, s) + (s — 53) 5(k — 1)} (17) 


(k+1)(k + 2)(k +3) 
Now applying the LT on boundary conditions we get, 


Hae dU(0,s)_ 1 4 dUQs) _ 4 
CO8)'= 0, dx s@+1 dx2 
= = 1 = 
Applying DTM, U(0,s) =0, U(1,s) = a4 and U(2,s)=0 
Put k = 0,1,2 ... in equation (17) we get, 
Ifk=0 20G3)=0 Ifk=1 *U(43)=0 
Ifk=2 « UG,8)=0 If k=3 -. U(6,s) =0 andsoon ... 
N = U k 
ow, U(x,s) >. U(x, s)(x) re 1% 
k=0 
Applying inverse Laplace Transform 
“ u(x,t) = xsint (Also shown in graphical presentation: see Fig.3) 
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6-30-20 m-20-10 m-10-0 m0-10 m10-20 20-30 
0-5 m5-10 10-15 15-20 20-25 = 25-30 


8-6 B64 O42 0-2-0 00-2 02-4 04-6 868 


Fig.1. Solution of u(x,t) = e*cost Fig.2. Solution of u(x,t) =x*+t  Fig.3. Solution of u(x,t) = xsint 


In this section, we have solved linear PIDEs with a convolution kernel of different order by using DLT, DET and 
LDTM. The graphical solution also been presented. It is observed that the DLT method reduces the linear PIDE 
with convolution kernel into algebraic equations (4), (9), and (14). Similarly, the DET method also reduces linear 
PIDEs into algebraic equations (5), (10) and (15). These algebraic equations take a lot of computational effort and 
time for their factorization and to achieve an exact solution of the problem. It has also been observed that the same 
solutions achieved by using LDTM method too. In the LDTM method, the Laplace transform reduces linear PIDEs 
into ordinary differential equations which are solved by the Differential Transform Method then the series solution 
of linear PIDEs is obtained by applying the inverse Laplace transform. In some cases, an exact solution can be 


achieved. 
5. Conclusion 


In this paper, the Laplace Differential Transform Method was successfully applied on linear Partial Integro 
Differential Equations with a convolution kernel. It is observed that, the solutions obtained by the LDTM method 
are the same as the DLT and DET methods. The LDTM gives a series solution and in some cases, the exact solution 
can be achieved. We firmly ascertained that LDTM takes less computational effort and time than others. It is 
concluded that LDTM is an alternate method to solve such kind of PIDEs. This method is a very simple, powerful 


and reliable for solving such problems. We hope some other types of PIDE can be used in various fields of 


modeling real life phenomena. 
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